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ABSTRACT 

The timescale for galaxies within merging dark matter halos to merge with each other is an important 
ingredient in galaxy formation models. Accurate estimates of merging timescales are required for 
predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar 
mass in central galaxies, and the statistical properties of satellite galaxies within dark matter halos. In 
this paper, we study the merging timescales of extended dark matter halos using iV-body simulations. 
We compare these results to standard estimates based on the Chandrasekhar theory of dynamical 
friction. We find that these standard predictions for merging timescales, which are often used in 
semi-analytic galaxy formation models, are systematically shorter than those found in simulations. 
The discrepancy is approximately a factor of 1.7 for M sat /Mhost ~ 0.1 and becomes larger for more 
disparate satellite-to-host mass ratios, reaching a factor of ~ 3.3 for M sa t/Mhost ~ 0.01. Based on 
our simulations, we propose a new, easily implementable fitting formula that accurately predicts the 
timescale for an extended satellite to sink from the virial radius of a host halo down to the halo's 
center for a wide range of M sa t/Mh os t and orbits. Including a central bulge in each galaxy changes 
the merging timescale by < 10%. To highlight one concrete application of our results, we show that 
merging timescales often used in the literature overestimate the growth of stellar mass by satellite 
accretion by rs 40%, with the extra mass gained in low mass ratio mergers. 
Subject headings: galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

As originally formulated by [Chandrasekhar (1943f), the 
deceleration of an orbiting point mass '"satellite," due to 
dynamical friction on a uniform background mass distri- 
bution, is given by 

4" "orb = -47TG 2 lll(A) M sat p h ost(< «orb) > 0) 

dt "orb 

where /O n ost(< "orb) is the density of background parti- 
cles with velocities less than the orbital velocity v rb of 
the satellite, M sat is the mass of the satellite, and A is 
the usual Co ulomb logarithm (e.g, Chandrasekhar] Il943l : 
lWhitdll97l . 

Accurate estimates of the effects of dynamical friction 
and the timescale for an orbiting satellite to lose its en- 
ergy and angular momentum to merge with a host are 
essential for many astrophysical problems. A thorough 
understanding of dynamical friction is particularly criti- 
cal for semi-analytic models (SAMs) of galaxy formation. 
The growth of dark matter halos via mergers can be com- 
puted either analytically or via dissipationless cosmolog- 
ical simulations, but the growth of galaxies depends on 
their dynamical evolution within larger dark matter ha- 
los. As a result, dynamical friction (and tidal stripping, 
which determines M sat ) provides a critical link between 
dark matter halo mergers and the galaxy mergers that 
determine, e.g., stellar masses, supermassive black hole 
masses, galaxy colors, and galaxy morphologies. Since it 
is computationally infeasible to simulate both dark mat- 
ter on cosmological scales and baryonic physics relevant 
to galaxies, SAMs will remain an important tool for in- 
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terpreting observations of galaxy formation and evolu- 
tion for the foreseeable future. 

There are, however, often significant uncertainties in 
directly applying equation ([1]) due to the approximate 
nature of the Chandrasekhar formula, and ambiguities 
in the appropriate value of the Coulomb logarithm and 
the definition of the satellite mass M sat for extended ob- 
jects (e.g., galaxies or dark matter halos). SAMs, for 
instance, typically use mild variations on the same basic 
formula for modeling the merger timescale r mcrgo induced 
by dynamical friction: 



Tdyn 



1.17 



/df Qorb M] 

In A Af„ 



host 



(2) 



(e.g., eq. 7.26 of iBinnev fc Tremaindll987t ). where or b 
contains information about the orbital energy and angu- 
lar momentum, /dt is an adjustable parameter, and Td yn 
is the dynamical timescale at the virial radius r v ; r of the 
host halo, related to the circular velocity at r v ; r , V^(r v i r ), 

by 



Vc( ' vir 

) \GM host 



1/2 



(3) 



Differences among different SAMs enter mainly in how 
each model treats In A, O or b, and /df, as well as how 
M sa t i s determined an d whe n r merge is determined 
(e.g.. iKauffmann et ail 119931: ISomerville fc PrimaclJ 
ll999tlCole et al.ll2000UCroton et al.ll2006h . Uncertainties 
in these parameters are reflected in equally large uncer- 
tainties in T morgc . 

The Coulomb logarithm, for example, should techni- 
cally be expressed as | ln(l + A 2 ); using In A is ap- 
propriate in the limit of large A (|Binnev k. Tremaind 
Il987f h The Coulomb logarithm represents the ratio of 
the largest and smallest impact parameters of field stars 
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that contribute to small-angle scatterings of the satel- 
lite: A = frmax/frmin- This motivates a conventional 
choice for the Coulomb logarithm: A = 1 + Mh st/Af sa t, 
where M sat is often taken to be the mass of the satel- 
lite when it enters the virial radius of the host halo 
of mass Mhost- Alternatively, the Coulomb logarithm 
can be taken to be a mass-in dependent constant (e.g., 
iCoorav fc Milosavlievid I2005D . When using equation 
(TTJ) to model the full orbits of satellites, another ap- 
proach is to assume that the Co ulomb logarithm varies 
over the course of the orbit (e.g., iHashimqto et al1l2003t 
IZentner et all 120051: iFellhauer fe Linl I2007D . These vari- 
ous prescriptions for the Coulomb logarithm can differ 
by a factor of 2-3. 

In addition, it is known that M sat cannot refer to just 
the bound mass for extended objects because tidally 
stripped material still in the vicinit y of the satellite 
also contributes to dy namical friction (|Fuiii et alj|2006t 
IFellhauer fe LinlfeOOTj ). The uncertainties in r me rge in- 
troduced by different prescriptions for In A and M sajt 
are important for galaxy formation models because 50% 
changes in the timescale for galaxies to merge within dark 
matter halos could significantly change the predictions 
for galaxy growth via merging or cannibalism, black hole 
merger rates, and the evolution of satellite galaxies in 
groups and clusters. 

Full dynamical models of the evolution of merging ha- 
los - complete with the physics of dynamical friction, 
tidal stripping, and gravitational shocki ng - have been 
creat e d by a number of groups fe.g.. iTavlor fc Babull 
20011 iTaffoni et ail l2003t IZentner et alj l2005t iGnedini 
20031) . These models, however, are based on ap- 
proximate treatments of the underlying physical pro- 
cesses, and obtaining merging timescales would often re- 
quire numerically integrat ing individual satellite orbits 
(| Velazquez fc White! 119991 ). In this paper we take a dif- 
ferent approach and instead compute merger timescales 
directly using numerical simulations. We focus on the 
range of satellite masses and orbital parameters that 
are relevant for dark matter halo mergers in hierarchi- 
cal galaxy formation models. Our goal is to provide a 
fitting formula that is as simple as possible but that still 
accurately reflects the physics of dynamical friction and 
tidal stripping. 

The remainder of this paper is organized as follows. 
Section 2 describes the numerical simulations we have 
performed. In Section 3, we investigate the merging 
timescales r mcrgc computed from our numerical simula- 
tions, propose a new fitting formula for r morge , and com- 
pare to previous results. Section 4 contains two sample 
applications of our proposed fitting formula: an estimate 
of the mass spectrum and orbital distribution of merging 
galaxies. Our results and conclusions, and their impli- 
cations, are reviewed in Section 5. Throughout this pa- 
per, we use "virial" quantities that are defined relative to 
200p c , where p c is the critical density. When necessary, 
we assume a cosmology of Q rn = 0.3, fl\ = 1 — VL m = 0.7, 
and H = 70 km s" 1 Mpc~ 1 . 

2. SATELLITE DECAY FROM NUMERICAL SIMULATIONS 

2.1. Description of simulations and initial conditions 

We have performed a number of numerical simulations 
to test the agreement between merger timescales pre- 
dicted by equation ([2]) and those derived directly from 



A^-bod y simulations. Ea ch of our simulations consists of 
a host iHernquistl (jl990f ) halo and a satellite Hernquist 
halo; the ratio of satellite to host mass and the initial 
orbital parameters of the satellite are varied from run 
to run. The host and satellite halos were constructed 
using AThost = 2 x 10 5 particles for the host and equal 
particle masses for the host and satellite particles, i.e., 
A^at = (M sat /Mhost) 2 x 10 5 . Each halo was constructed 
by sampling the full phase-space distribution function 
(under the assumptions of spherical symmetry and veloc- 
ity isotropy) and was tested to be stable when evolved in 
isolation. In addition, two simulations included a (self- 
consistently constructed) Hernquist bulge in both the 
host and satellite to test the effects of baryonic com- 
ponents on merging timescales. 

The numerical simulations were run using GADGET- 
2 (lSpringe]|l2005h . Outputs were saved every 0.1 Gyr, 
or equivalently, sa 0.06 Td yn . To ensure that our simu- 
lations are not affected by numerical artifacts, we have 
performed several convergence runs with respect to both 
the particle number (up to 10 6 particles in the host halo) 
and the force softening e (which is reduced by up to a fac- 
tor of 5 from our fiducial value of e = 10~ 3 r v ; r ). Because 
our simulations involve only gravity, the length, mass, 
and time scales can be rescaled for any host halo mass; 
for definiteness, however, we quote our results for host 
halos of total mass A/host = 10 12 M q and Hernquist scale 
length a = 40 kpc. Matching to a iNavarro et all (|1997l 
hereafter NFW) profile with a virial mass equal to the 
Hernquist halo's total mass and an identical density at 
r -C a, our standard host halo correspon ds to an NFW 
halo w ith a concentration « 8.5 (see, e.g., Springel et al. 
I2005al for more details on how to relate Hernquist and 
NFW halos). A comparison run using an NFW halo had 
a merging timescale within 5% of that from the corre- 
sponding run with a Hernquist halo. 

In the absence of dynamical friction, the orbit of a test 
particle in a static spherical potential is entirely specified 
by its energy E and angular momentum J. An equivalent 
parametrization for bound orbits is to use the orbital 
circularity 

V = ~~7T^ ' (4) 

Jc{E) 

which is the specific angular momentum relative to the 
specific angular momentum of a circular orbit with the 
same energy [rj is related to eccentricity e by rj — (1 — 
e 2 ) 1 / 2 ], and r c (E), which is the radius of a circular or- 
bit with the same energy as the particle in question. 
Since dynamical friction dissipates energy, the orbit of 
a satellite galaxy in a host galaxy generally depends on 
its initial position r(t = 0) in addition to E(t = 0) and 
J(t = 0). 

We explore a range of orbital circularity n, energy 
r c (E), and satellite-to- host mass ratios M sat /Mh os t in 
our simulations. A summary of the production runs and 
the resulting r me rge is presented in Table Q] The ranges 
of the parameters are chosen so that the satellites can 
plausibly undergo significant orbital evolution within a 
Hubble time (see Sec. 13.21 for more details). Systems 
with mass ratios M sat /Mh os t much smaller than 1:40 are 
therefore dynamically uninteresting and excluded. 

2.2. Defining timescales: angular momentum loss 
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TABLE 1 
List of Simulations 



Run M sa t/M host circularity r] r c (E)/r 

vir Emerge 



(1) 


(2) 


(3) 


(4) 


(5) 


3b 


0.3 


0.78 


1.0 


4.4 


3d 


0.3 


1.00 


1.0 


6.9 


5b 


0.2 


0.46 


1.0 


3.5 


5c 


0.2 


0.65 


1.0 


4.5 


5d 


0.2 


0.78 


1.0 


6.0 


5e 


0.2 


1.00 


1.0 


9.5 


10a 


0.1 


0.33 


1.0 


4.85 


10a* 


0.1 


0.33 


1.0 


1.3 


10b 


0.1 


0.46 


1.0 


6.0 


10c 


0.1 


0.65 


1.0 


8.25 


10c* 


0.1 


0.65 


1.0 


7.8 


lOd 


0.1 


0.78 


1.0 


10.75 


14b 


0.0707 


0.46 


1.0 


9.0 


14c 


0.0707 


0.65 


1.0 


12.5 


20a 


0.05 


0.33 


1.0 


10.0 


20b 


0.05 


0.46 


1.0 


11.0 


20c 


0.05 


0.65 


1.0 


17.0 


40b 


0.025 


0.46 


1.0 


22.0 


40c 


0.025 


0.65 


1.0 


38.0 


lOal 


0.1 


0.33 


0.65 


3.5 


10a4 


0.1 


0.33 


0.80 


1.3 


lObl 


0.1 


0.46 


0.65 


1.3 


10b4 


0.1 


0.46 


0.80 


5.0 


lOcl 


0.1 


0.65 


0.65 


5.6 


10c2 


0.1 


0.65 


0.70 


5.5 


10c3 


0.1 


0.65 


0.75 


6.2 


10c4 


0.1 


0.65 


0.80 


6.35 


10c5 


0.1 


0.65 


0.90 


7.5 


Note. 


— Description 


of columns: 







(1) Name of simulation 

(2) Ratio of initial satellite mass to initial host halo mass 

(3) Initial orbital angular momentum, parametrized by circular- 
ity 7) 

(4) Initial orbital energy, parametrized by r c (E)/r v - lr 

(5) Dynamical friction merging time r me rge, in Gyr, for a host 
with virial mass M^ost = 10 12 Mq, measured from the simula- 
tion as described in Sec. [2] 

(* indicates run with baryonic component; see ^12.31 for details) 

The timescale for a satellite to merge with its host halo 
can be defined in a number of different ways. One com- 
mon definition is to take a fiducial radius of the bary- 
onic component assumed to reside at the center of the 
halo and assume that the satellite has merged when its 
separation from the host's center equals this radius. In 
this work, we instead consider a satellite merged when 
it has lost all of its specific angular momentum j = rvt 
(relative to the host). In practice, this definition agrees 
with the more commonly-used one in most cases but also 
works well for situations in which the definition based on 
orbital separation can give undesired results. Consider 
highly eccentric orbits, in which the satellite can come 
very close to the center of the host while retaining sig- 
nificant orbital energy; satellites on these orbits often do 
not merge for multiple dynamical times following the first 
close encounter. The specific angular momentum of the 
satellite should, however, be a non-increasing function of 
time, so using j(t)/jo [where jo = j(t = 0)] is therefore 
our preferred definition for merging. 

Fig. Q] shows the trajectories (top) and angular mo- 
mentum decay (bottom) for three sets of live-satellite 
simulations starting from r v i r , with M sat /A/host = 
0.025, 0.05, and 0.10, for two different orbital circular- 
ities, 77 = 0.46 (solid curves) and 77 = 0.65 (dashed 
curves). For a given circularity, systems with more dis- 
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Fig. 1. — Top: Trajectories of satellites (i.e., separation between 
satellite and host centers) from six simulations with two different 
orbits — T] = 0.46 (solid) and 0.65 (dashed) - and three different 
mass ratios: 1:10 (black), 1:20 (blue), and 1:40 (green). For all 
cases, satellites on more eccentric orbits (lower r)) or with higher 
mass ratios merge more quickly. Bottom: angular momentum de- 
cay for the same satellites. The loss of angular momentum ini- 
tially correlates well with pericentric passages but at later times 
dj/dt ~ constant. 

parate masses merge much more slowly, as is expected 
from basic dynamical friction considerations. For a given 
mass ratio, the satellite on the more eccentric orbit (i.e. 
smaller 77) loses angular momentum faster, resulting in a 
more rapid merger. The difference is non-negligible even 
for mass ratios as similar as 1:3. In this satellite 
with 77 = 0.78 merges in approximately 4.4 Gyr while a 
satellite on a circular orbit takes over 50% longer - 6.9 
Gyr - to merge (see last column of Table 1). 

Fig. [T] highlights the angular momentum loss process. 
A comparison between the top and bottom panels shows 
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that the bulk of the angular momentum loss initially co- 
incides with pericentric passages (e.g., at 1 and 4.5 Gyr 
for the 1:20 blue curves) and t he accompanying tidal 
stripp ing and shocking (see also iBovlan-Kolchin fc Mai 
20071). It is also interesting to note that while the an- 
gular momentum loss is initially somewhat impulsive, it 
later (t > 5 Gyr) becomes nearly constant in time. Phys- 
ically, the approximate constancy of dj/dt ~ j/T merg e(r) 
follows from the fact that the merger timescale at any 
radius is roughly linearly proportional to radius (see eq. 
[5] discussed below), as is j itself. It is also interesting 
to note that a satellite can continue to orbit for several 
giga-years even after its initial orbit has dipped within 
5% of r vil of the host (e.g., the 1:20 merger in Fig. [T]). 
In fact, at the first pericentric pass within ~ 5% of r v i r 
(which occurs at t w 6.5 Gyr for rj = 0.46), the satellite 
has only lost 60% of its initial angular momentum. Such 
satellites should not be considered merged even though 
they may pass very close to the center of the halo (within 
the radius at which a central galaxy might lie); our cal- 
culations show that these satellites will orbit for a signif- 
icantly longer time before actually merging. 

When the mass of the satellite is much less than the 
mass of the host, r mer ge becomes prohibitively long to 
study with simulations, as numerical relaxation becomes 
an important effect. For these cases (only three runs: 
20c, 40b, and 40c), we determine r morgG by linearly ex- 
trapolating j(t)/jo. This tends to be a reasonable mea- 
sure of the merging time as long as we extrapolate af- 
ter the first two pericentric passages, where gravitational 
shocks and tidal stripping cause substantial changes in 
the angular momentum. After the first two pericentric 
passages, the angular momentum loss occurs at a roughly 
uniform rate. 

2.3. Effects of baryons 

Although the ratio of a galaxy's stellar mass to the 
virial mass of its dark matter halo does not exceed w 0.1, 
it is conceivable that neglecting baryons could lead to 
significant errors in the predicted merging timescales. In 
particular, bulges in elliptical galaxies are significantly 
more dense than dark matter halos at the same physical 
scale, meaning that bulges are more resistant to disrup- 
tion via tidal shocking. Thus, the expectation is that 
merger timescales would be shorter for full galaxy mod- 
els than for dark matter halos alone. 

In order to test the effects of the baryonic components 
of galaxies on merging timescales, we ran two additional 
simulations. Both simulations used M sat /M host = 0.1, 
but we included stellar bulges with M+/M dm = 0.05 in 
both the host and satellite for each simulation. 4 The 
effective radii of the bulges were chosen to be represen- 
tative of those measured bvlShen et al.l (120031 ) from the 
Sloan Digital Sky Survey (lYork et al.ll2000ri . JL cx M° 56 , 
and the orbits of the two runs are identical to runs 10a 
and 10c. 

In Fig. [2j we compare the angular momentum loss as a 
function of time for the runs without stars (solid curves) 
and with stars (dashed curves). Including stars does lead 

4 The bulge + dark matter halo systems for the initial conditions 
are set up self-consistently and are ve ry stable over many dynamica l 
times when evolved in isolation; see IBovlan-Kolc hin & Ma (2007) 
for details. 
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Fig. 2. — Effects of including baryonic components. The solid 
curves show the angular momentum loss for runs 10a (left; r\ = 
0.33) and 10c (right; r\ = 0.65). The dashed curves show j(t)/jo 
for runs that have a stellar bulge (with M v i T / M^ u i ge = 20) in both 
the host and satellite but are otherwise identical. Including baryons 
leads to merger timescales that are < 10% shorter than those of 
the corresponding dark matter-only simulations. 

to faster mergers, as expected, but the difference in merg- 
ing timescales is less than 10%, which is quite modest. 
When the satellite first enters the host's virial radius, the 
bulges contribute negligibly to dynamical friction drag 
because the baryons contribute < 10% to the system's 
total mass. At late stages in the merger, the bulges may 
be more important, as the dark matter halo is more ef- 
ficiently stripped than the bulge. The satellites spend 
the vast majority of their time at large radii, however, 
so that the error introduced by ignoring baryons is small 
when considering the total time required to merge from 
the virial radius (see Fig. [5]). 

3. COMPUTING MERGER TIMESCALES 

In this section we provide a simple parameterization of 
the merging timescales determined from the numerical 
simulations discussed above. Comparisons between this 
new fitting formula and some commonly used ones in the 
literature are also presented (Sec 3.3). 

3.1. A fitting formula 
We fit the merger timescales from the simulations in 



Table 1 to a simple formula, 



' merge 
7"dyn 



A 



(M host /M sat ) b 
ln(l + Mhost/M sa t) 



exp 



ME) 



( 5 ) 

where the constants 6, c, and d parameterize the depen- 
dence of the merger timescales on the host-to-satcllitc 
mass ratio Mh os t /M sat , orbital circularity n — j/j c (E), 
and orbital energy r c (E), respectively. Both n and r c (E) 
are computed self-consistently using the Hernquist po- 
tential, not using the two-body approximation. To avoid 
ambiguities in the definition of satellite and host mass 
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as a function of the location of the satellite (or time), 
we define the masses Mh os t and M sa t here as the virial 
masses of the host and the satellite when entering the 
host's virial radius. For the same reason, we use the spe- 
cific angular momentum j rather than the full angular 
momentum J in computing the orbital properties of the 
satellites; since the circularity depends only on a ratio of 
angular momenta, this choice does not affect calculations 
of 77. The dynamical time Td yn is the dynamical time at 
the host's virial radius r v ; r given by equation ([3]). Note 
that Tdyn = 0.1 -ff -1 (where H is the Hubble constant) 
independent of the host halo mass; at z = 0, Td yn = 1.4 
Gyr for H = 70kms~ 1 Mpc -1 . 
With the definitions above, we find that 



A = 0.216. 



1.3, 



1.9, d=1.0 



(6) 



provide fits to the simulation results with a standard de- 
viation of < 7% and a maximum deviation of 12.5%, 
a significant improvement over equation ([2]). To ob- 
tain the values in equation (|6|), we have chosen to fit 
to the parameter d separately from A, b, and c. For 
A, b, and c, we have fitted to r mcrgo from the subset 
of simulations in which the satellites begin at the virial 
radius of their host with an orbital energy E — E c (r v i r ), 
i.e., with r c (E) = r v i r , and the fit is independent of 
d. To determine the value of d in equation ([6]), we 
use the series of simulations with -Msat /Mhost = 0.1 and 
r c (E)/r vil = [0.65, 0.7, 0.75, 0.8, 0.9, 1.0] (all starting at 
r vir ). Holding A, b 7 and c fixed, we investigate whether 
any single choice of d results in a consistently accurate 
prediction of the merging timescale. 

The close agreement between the numerical simula- 
tions and the fits in equation ([5]) is illustrated by the cir- 
cles in Fig. [3] and squares in Fig. The two figures show 
the ratio of the timescale measured from each simulation, 
T merg e(sim), to the timescale predicted by equation ([3]), 
T"merge(fit), for a wide range of orbital energy (0.65 < 
r c (E)/r V i r < 1.0), orbital circularity (0.33 < 77 < 1), and 
satellite-to-host halo mass (0.025 < M sat /M\ lost < 0.3). 
In addition, Fig. [3] shows that using d — 1.0 (circle sym- 
bols) in equation ([5]) results in good agreement between 
predicted and measured r mcrgo , while the often used 
d = 2.0 (diamond symbols) results in a poor match, sys- 
tematically underestimating T mcr gc (see §3.3.21 for more 
discussion). In Fig. [H we also compare T merg e(fit) to 
that predicted by a fiducial SAM, r morgG (SAM); this is 
discussed in ^3.31 below. 

A further point of interest is that using da; 1 in Equa- 
tion ([5]) also provides a significantly better fit than d = 2 
when considering the merging timescale as a function of 
position along a given orbit. This suggests that equa- 
tion ([5J can be used for starting radii other than r v - u , 
though its accuracy will certainly decline for r <C r v ; r . 
This possibility is not explored in any more detail in this 
paper, however, and we restrict ourselves to starting radii 
of r vir . 

3.2. Range of validity 

We have chosen to simulate parameters that corre- 
spond to probable satellite orbits in a ACDM cosmol- 
ogy. The mass ratios we have considered - 0.025 < 
M sat /Mhost < 0.3 - cover a reasonable range, as lower 
mass ratios are unlikely to merge within even multiple 
Hubble times, while mass ratios near unity merge on 
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Fig. 3. — Comparison of merger time measured in our iV-body 
simulations, r mcrgc (sim), with the predicted T mC r ge (nt) based on 
eq. J5J, as a function of orbital energy r c (E)/r v i r . The orbital 
circularity is color-coded: n = 0.33 (red), 0.46 (green), and 0.65 
(black). The circles use the best- fit parameters in eq. ((Sj and show 



the close agreement in t u 



between the simulations and the fits. 



The diamonds, on the other hand, illustrate the poor match when 
the energy dependence in eq. J5]l is taken to be the conventional 
d = 2.0 instead of d = 1.0. All points are for M sa t/Mh 03 t = 0.1. 

nearly a dynamical time. We note that Eqn. [5] does in- 



deed tend toward r„ 



Tdyn for M sat = Mhost and 



77 ss 0.5; it therefore should give a reasonable approxima- 
tion even for our untested regime of M sa t/Mhost > 0.3. 

We have covered a range of circularities - 0.3 < r\ < 1.0 
- that includes the most likely values based on analy- 
ses of orbits i n dark matter simulations (77 ~ 0.5; e.g., 
Benson! [2fJ05t IZentner et ail 120051 : iKhochfar fc Burkertl 
2006h . Equation ((5]) should not be used for rj < 0.2 
because such orbits pass sufficiently close to the center 
of the host halo on their first pericentric passage that 
interaction with the central galaxy is important and will 
decrease the merging timescale. 

The range of orbital energies considered here - 0.65 < 
r c (E)/r v i r < 1.0 - also covers the peak values of distribu- 
tions seen in cosmological simulations. We have limited 
ourselves to one specific scaling of dark matter halo scale 
radii - a cx M 5 - that is motivate d by analyses of cos - 
mological TV-body simulations (e.g.. lBullock et aLll2001h . 
which show t hat c oncentrations scale as c cx M~ - 13 . 
Taffo ni et al.l {2003) found that the dependence of r mcrgc 
on the satellite concentration is relatively weak: using 
their eq. 27, varying c sa t/chost between 1 and 2 results in 
a change in r mcrgo of only « 20%. 

3.3. Comparison to Previous Work 
3.3.1. Exponent b: dependence on mass ratio 

One important difference between equations ((2]) and 
([5]) is the non-linear dependence on the mass ratio in 
equation ((5]). Our fitted value of b = 1.3 reflects an 
important difference between point masses and realis- 
tic dark matter satellites sinking in a larger host po- 
tential, as the dark matter halos continually lose mass 
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Fig. 4. — Comparison of merger times computed from Af-body 
simulations, fitting formula eq. J5J, and a fiducial SAM using eq. 
and two choices of Coulomb logarithm (black and blue curves) , for 
a variety of satellitc-to-host mass ratios M Ba ,t/M^ oat and orbital 
circularities n, with r c (E) /r v i r = 1. The square symbols show 
the ratio T me rge(sim)/T me rgc(fit) and illustrate that the simulations 
and fit agree well, with an error of < 10%. The orbital circularities 
are color-coded: n = 0.33 (red), 0.46 (green), 0.65 (black), 0.78 
(magenta), and 1.0 (yellow). The fiducial SAM, on the other hand, 
generally underestimates the merging time, particularly at small 
M sa t/Mh oa t, as shown by the three curves for different circularities: 
r) = 0.25 (solid), 0.5 (dashed), and 0.75 (dot-dashed). The SAM 
here assumes f^f = 1, O rb = n°' 7s , and r c (E)/r v i r = 1. The 
different color curves correspond to different choices of Coulomb 
logarithms for the SAM: ln(l + M^ ost / 'M sa t) [black curves] and 
| ln(l + M h 2 ost /A/ s 2 at ) [blue curves] 

via tidal stripping and gravitational shocking against the 
host potential. The effective mass of the satellite over 
its lifetime is therefore smaller (sometimes substantially- 
smaller) than its mass when entering the host halo. It 
is therefore not surprising that extended satellites sink 
more slowly than the standard formula predicts, leading 
to b > 1. 

The exact exponent b turns out to be difficult to model 
from first princi ples as it depends o n a n umber of factors. 
As an example, iFuiii et aLl (|2006D and iFellhauer fe Lin! 
(|2007f ) have shown that the effective mass (with regards 
to dynamical friction) of a sinking satellite is actually 
larger than the instantaneous bound mass because the 
mass that has been recently stripped also contributes to 
the drag force on the satellite. We see the same effect in 
our simulations, significantly complicating any effort to 
assume that either the bound mass or the instantaneous 
tidal mass is the relevant mass at every time. 

3.3.2. Parameters c and d: dependence on satellite orbits 



A common 
dependence 



choice for 



(S 



m equat ion 



tiack 
3jTfoi 



the 

M- 



angular 



is 



Somerville fc Primact : 1999 nCole~ et al 



momentum 
oc T7 - 78 
which 



2000), 



Lacev fc Cold (|1993|) found to be a good match to their 
integration of the orbit-averaged equations for energy 
and angular momentum loss for a point mass subhalo 
due to dynamical friction in an isothermal potential. 



Ivan den Bosch et alj (1999) found a somewhat weaker 
angular momentum dependence in their numerical sim- 
ulations (using point-mass sat ellites and an isothe rmal 
potential), r mergc oc 77 053 , while [Taffoni et all (|2003f ) pro- 
vide a more complex fitting formula for both point-mass 
and live satellites with realistic internal structure in an 
NFW potential. 

We find that an exponential dependence - r mcrgo oc 
exp(c?7) with c = 1.9 - provides a better fit to our sim- 
ulations than a pure power law: a fit in which we force 
TVnerge oc ry - 78 has a standard deviation of 14 % and a 
maximum deviation of 42% when compared to the simu- 
lations (compared to 6.7 % and 12 %, respectively, when 
using our fit). It is important to note, however, that it is 
the long Tmerge for the 77 = 0.78 and 77 = 1.0 simulations 
that result in better agreement using the exponential fit 
rather than the power-law fit. If we restrict our attention 
to runs with 0.33 < 77 < 0.65, we find that r mcrgc tx r/ - 7& 
does indeed provide a good fit to our simulation results. 

Our revised fitting formula has a somewhat weaker de- 
pendence on 77 for 77 < 0.5 and a stronger dependence 
for 77 > 0.5 relative to that found by previous work with 
point-mass satellites. Both of these limits are plausible: 
for small 77 (nearly radial orbits), a slight change in 77 
should have a negligible effect on r mcrgo since the orbits 
will have small pericentric distances. For large 77 (nearly 
circular orbits), live satellites still lose mass at large radii, 
reducing their mass and lengthening their r mcrgc relative 
to an equivalent point-mass satellite. The same is true 
for live satellites on more radial orbits but to a lesser 
degree, as angular momentum losses correlate well with 
pericentric passages (see Fig. [T]); the net effect creates a 
stronger dependence on 77 for 77 > 0.6 for live satellites 
than for point-mass satellites. 

For the orbital energy dependence of the merger 
timescale, we find that d = 1.0 is the best fit to 
our simulations. This is to be contrasted with the 
canonical value of d 2.0 in prior work. The lat- 
ter ^is_j 1 2p_rop_riate_JcT anjsothermal host (e.g., eq. 7.26 
of iBinnev fc Tremainel Il987f ). or for point mass satel- 
lites sinking on circular orbits in a static potential, as 
can be determined by integrating equation |T]). This 
yields T merg e < x (rc (E)/r vil -) 1 - 97 for an NFW potential 
(Ta ffoni et all I2003f) . We find that for point -mass satel- 
lites on circular orbits in a host halo with a Hernquist 
density profile and scale radius a, the energy dependence 
of T merge has an analytic expression, 



' merge 



1 

To 



1 

3 X 



1 



tan 



(7) 



where x = r c (E)/a. Equation ([7]) can be well- 
approximated by T merg0 oc r c (E) 2 - 3 for 0.2 < x < 30, 
comparable to d = 2. The difference between these re- 
sults for point masses and our best-fit value of d = 1 
lies in our study of live satellites and in our definition of 
Mgat. Because of ambiguities in defining M sat at small 
radii when significant tidal stripping has occurred, we 
have chosen to use virial quantities exclusively in equa- 
tion (O, including for the satellite mass M sat . Were we to 
use Zoca^7/-defined quantities such as M sat (£), we would 
expect that M sa _ t (t) oc r(t) for an isothermal sphere, so 
that the merger time evaluated as a function of radius 
is given by r mcrgo oc r 2 (t) / M sat {t) oc r(t). It is therefore 
not surprising that we find d ~ 1 to be a better match 
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Fig. 5. — Contours of equal r mcrgc in the space of satcllitc- 
to-host mass ratio Msat/Mhost an d orbital circularity r], assuming 
r c (E) = r v i r . The black solid curves represent our fitting formula 
eqs. J5j and J5J; the blue dashed curves correspond to the same 
SAM shown in Fig. [4] The labelled times are for a host halo with 
Mhost = 10 12 Mq and a = 40 kpc. The SAM generally under- 
predicts r mcrgc , with the largest discrepancy at large r\ and small 
M sat /M host . 

than d 2 to our numerical results. 

3.3.3. Numerical comparison 

In Figs. 2] and [5] we quantify the differences between 
equation ([5]) and prior work discussed in the previous two 
subsections. Fig. [4] compares r mcrge (SAM) from a fidu- 
cial SAM using equation {SJ and T mcrgc (fit) from our fits 
to equation ([5]) for a range of satellite-to- host mass ratio 
Msat/Mhost and orbital circularity (f) = 0.25, 0.5, and 
0.75 for the solid, dashed and dot-dashed curves, respec- 
tively). The SAM here assumes f df = 1, O orb = ry 78 , 
and In A = ln(l + Afhost/M sat ). The curves illustrate that 
the fiducial SAM underestimates the merger timescale at 
all masses we have considered regardless of the orbital cir- 
cularity. Moreover, the disagreement grows as the satel- 
lite mass decreases relative to the host: at Msat/Mhost = 
0.2 the SAM timescale is too small by approximately 
50%, while the discrepancy increases to a factor of 3 at 
Msat/Mhost = 0.025. If a Coulomb logarithm of A = 
1 + (M hos t/M sa t) 2 is used (|Somerville fe Primacklll999f ). 
the discrepancy increases further by approximately a fac- 
tor of 2 for M sat <C Mhost- In order to provide a better 
match to simulations or observations, some groups ad- 
just the normalizat i on of T mer ge by s etting fdf ^ 1 (e.g. , 
Ide Lucia fe Blaizotl [20071 Ut = 2; iBower et al.l [20061 
= 1.5; and iNaeashimaet all [20051 fdf = 0.8). Al- 
though fdf > 1 goes in the right direction, Figure [4] shows 
that a constant fdf does a poor job of matching the nu- 
merical results. 

Fig. compares the contours of constant r morgo in the 
space of mass ratio M sat /Mh os t and orbital circularity rj, 
as predicted from our fits in equation ([5]) (solid curves) 
and the same SAM from Fig. 0] (dashed curves). The 
two sets of curves are for the same five merger times: 2, 



5, 10, 15, and 20 Gyr. All orbits that lie to the right 
of and below a given contour will merge within the time 
corresponding to that contour, while those to the left of 
and above a contour will not. Fig. [5] shows that even 
the predictions for orbits with the fastest merging times 

- those with low circularities and large mass ratios (the 
lower-right portion of the figure) - differ between the 
SAM and our fit. At larger circularities or smaller mass 
ratios, the two predictions deviate significantly. As an 
example, the prediction from our revised formula for 20 
Gyr merging times is similar to the curve for 10 Gyr 
merging times in the SAM. 

Our results can be compared to t he fitting for- 
mu lae provided by iTaffoni et alj ((2003) [and updated 
by iMonaco et all 12007} . which are valid for 0.3 < 
r c (E)/r v i T < 0.9. The overall trend we find is qualita- 
tively similar to the results of Taffoni et al. (e.g., their 
Fig. 7), but the quantitative details are somewhat differ- 
ent. We find that massive satellites (M S at/Mh os t ^ 0.05) 
sink more slowly than predicted by Taffoni et al. while 
light satellites sink more quickly. For example, we pre- 
dict a timescale that is 25% shorter than Taffoni et al. 
for Ms a t/Mh st = 0.025, r c (E)/r vil = 0.75, and r) = 0.65 
but find a timescale that is 40% longer for a satellite 
with M sat /Afhost = 0.1 and identical orbital parame- 
ters. Fig. 7 of Taffoni et al. also shows that the de- 
viation of live satellite merger times from a fiducial SAM 
prediction is not monotonic with M sat /Mh os t as is sug- 
gested by our Fig. [4] and equation ((5|). In particular, 
for Msat/Mhost < 0.001 and M sat /M host > 0.2, 

Emerge 

is approximately the same in the extended and point 
mass cases. We have not explored mass ratios lower than 
1:100, however, because the sinking times exceed many 
Hubble times. 

After this paper was submitted, iJiang "eTall ([2007L 
hereafter J07) independently investigated the merging 
timescales of satellites using A-body simulations and 
came to a similar overall conclusion as that in this pa- 
per, namely, that Eqn. [2] generally underpredicts merg- 
ing timescales. J07 used a different .A-body approach 

- cosmological simulations including hydrodynamics - 
and obtained a fitting formula that differs in its de- 
tails from ours. Specifically, J07 found b — 1 (rather 
than 1.3) and a significantly weaker dependence on cir- 
cularity. The aforementioned differences between the 
methods used in this study and J07, as well as differ- 
ences in definitions (e.g., we define halo radii relative to 
200 p c while J07 use A v (z) p c ) make a direct comparison 
non-trivial. However, even for our most major mergers 
(Msat/Afhost > 0.1, runs 3x, 5x, and lOx), we still find a 
significantly stronger dependence on r\ than J07: whereas 
the fitting formula in J07 predicts a difference in r merge 
of only 30% between rj = 0.46 and rj = 1 orbits, we find a 
difference of 160% (runs 5b and 5e). In addition, we find 
that at fixed orbital parameters, the difference in merger 
times between our most major merger simulations (the 
0.3:1 and 1:5 runs) is somewhat better fit by b = 1 than 
b = 1.3. Most of the non-linear dependence on mass-ratio 
that we find in equation 6 is thus due to the larger mass 
ratio (more minor) mergers. This suggests that the dif- 
ference between our preferred value of b = 1.3 and J07's 
value of b = 1 may be due to the range of mass ratios 
considered. 
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Fig. 6. — Fraction of stellar mass gained through mergers (i.e., 
satellite accretion) at the host's virial radius r v ; r that makes it 
down to the central galaxy in time t, as a function of the satcllite-to- 
host stellar mass ratio . The thick solid line is the number of halo- 
halo mergers at r v i r as a function of mass ratio (normalized to 
integrate to unity). We consider the predictions of both the fiducial 
SAM from Fig.[I]and our eqn. J5]l. The fiducial SAM predicts that 
approximately 40% more mass is accreted in one Hubble time, all 
of which is in low-£* mergers. 

4. A SAMPLE APPLICATION 

The merger timescale given by equation can be 
applied to a wide variety of astrophysical problems. We 
defer most of these to future papers. To give one concrete 
example, however, we briefly consider the dissipationlcss 
growth of stellar mass for galaxies at the centers of dark 
matter halos. 

In a finite time period, whether or not galaxies can 
merge following the merger of their much more extended 
host halos depends on r mer ge due to dynamical friction. 
Given the dark matter halo-halo merger rate d 2 N m /d(,dz 
at some redshift z, defined as the average number of 
mergers with mass ratio £ = M2/M1 < 1 leading to 
halos of mass M — M\ + M2 per halo per redshift per 
£, the corresponding growth of stellar mass over a time 
interval t can be computed using 



d 4 M* 



M*(M 2 ) 



d 2 N m d 2 P 



e(t 



(8) 



d£dr]dEdz v ' d£dz drjdE 

where d 2 P/drjdE is the probability of two dark mat- 
ter halos merging with orbital circularity 7/ and energy 
E (typically measured when the two halos' virial radii 
touch), is the unit step function, and M*(M2) is the 
stellar mass as a function of dark matter mass. The step 
function in equation §8$ enforces the fact that only galax- 
ies in halos with r mer ge < t can merge with the central 
host galaxy in time t. A full calculation of dM+/dt would 
require an appropriate integral of equation ([5]) over time 
t, so that all halo-halo mergers are accounted for. Re- 
sults of the full calculation will be presented in Boylan- 
Kolchin et al. (in preparation) ; here we consider instead 
the simpler problem of using equation ([8]) to calculate 



the fraction of the stellar mass accreted at redshift z at 
r v ir that makes it down to the center within time t. In 
particular, by integrating equation 1(5)), the stellar mass 
accreted in satellites of different mass, circularity, or en- 
ergy can be computed. 

We first consider the mass spectrum of merging ob- 
jects that determines whether most of the stellar mass 
added to the galaxy comes from a small number of mas- 
sive progenitors or a large number of low-mass progen- 
itors. This is obtained from equation {5| by integrat- 
ing over 77 and E, taking into account the dependence of 
Tmerge on each of these orbital parameters. IZentner et al.1 
(2005) found that dP/d-q oc i] 12 (1 - r)) 1 - 2 , which we 
use in equation {5}. For simplicity, we assume that all 
orbits have an energy equal to the energy of a circu- 
lar orbit at the host's virial radius (r c (E) = r v i r ), i.e., 
dP/dE = S[E — E c (r v - U )]. The halo merger rate is taken 
from the analysis of Fakhouri & Ma (2007, in prepara- 
tion), who have analyzed merger rates in the Millennium 
simulation (|Springel et al. 2005b) and find that 



d 2 N m 
d£dz 



-1.92 



exp 



0.125 



0.4 



(9) 



Finally, we assume that dark matter halo mass and 
the stellar mass of its central galaxy scale as Md m oc 
AfJ- 5 , which i s appropriate for mass ive galaxies in groups 
and clusters (iGuzik & Seliakl l2002t iHoekstra et~alj |2005l; 
iMandelbaum et al.ll2006h . 

Fig. [6] shows how different dynamical friction prescrip- 
tions lead to different cutoffs at the low mass end for 
the fraction of mass assembled via mergers of a given 
(stellar) mass ratio £*. We plot both the predictions of 
our numerical fit in equation (|5|) and the fiducial SAM 
from Fig. 3J Techincally, we should consider orbits with 
r] < 0.2 separately from those with 7/ > 0.2, as the dy- 
namics of very radial mergers can be strongly affected by 
the stellar components of the merging galaxies at first 
pericentric passage, introducing complications that we 
have not considered in our current analysis. For illustra- 
tive purposes, however, these orbits, which account for 
ss 9% of the accreted mass at r v i r , are included in Fig. [BJ 
In all of the calculations in Fig. [6l a wide range of 
~ 0.1 — 1 contribute to the mass growth of the central 
galaxy. For a canonical SAM dynamical friction formula, 
however, a non-negligible amount of extra mass is added 
in low mass ratio mergers, relative to the predictions of 
our formula for r mcrg0 . In fact, dM+/dt is approximately 
40% larger for t = tu for the canonical SAM formula 
for T mcrgc , and all of the extra mass is accreted in low-£ 
mergers. 

Dynamical friction also distorts the orbital distribution 
of satellites at r v ; r of the host, preferentially selecting low 
angular momentum orbits for accretion onto the central 
galaxy. This effect is shown in Fig. [71 which shows the 
initial distribution of orbital circularities (dashed curve; 
taken from Zentner et al. and assumed to be inde- 
pendent of ^) along with the distribution of circulari- 
ties for satellites merging within a given fraction of a 
Hubble time (r mcrgc /7j ff = 0.1, 0.316, 1.0, and3.16). For 
Tmerge/ tu 3> 1, the stellar mass accreted as a function of 
circularity should approach the circularity distribution at 
r v i r , as all satellites will eventually merge. As T mcTge /tH 
decreases, however, the peak of the distribution shifts to 
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Fig. 7. — Impact of dynamical friction on the orbital properties 
of satellites that merge with a central galaxy. The input distribu- 
tion of circularities r/ for halo-halo mergers at r v i r with $ > 10~ 3 
is shown with the dashed curve (dP/dr], taken from Zent ner et al.l 
120051) . The solid curves show the predicted distribution of r\ for 
satellite galaxies that merge onto the central galaxy within a spe- 
cific fraction of a Hubble time using equation J5}. Higher mass 
ratio and lower angular momentum orbits preferentially accrete, 
leading to a shift in the dP/drj distribution to lower circularities. 

smaller and smaller circularities because low angular mo- 
mentum satellites merge more rapidly. Fig. [7] shows that 
the peak of the input distribution of circularities (dP/drj) 
is at rj — 0.5 but that most of the mass comes in via satel- 
lites on orbits with rj = 0.42 at r mer ge/ii/ = 1 and with 
?7 = 0.30 for T mclgc /tH = 0.1. Shifting the peak of the 
circularity distribution from 0.5 to 0.42 corresponds to 
a 30% reduction in the pericentric distance of the or- 
bit, which can have a significant e ffect on the properties 
of th e merger remnant (see, e.g., IBovlan-Kolchin et al.l 
120011 . 

5. CONCLUSIONS 

Applications of dynamical friction to problems in 
galaxy formation are often based on classic results de- 
rived for point masses sinking in a galaxy with an isother- 
mal density profile (e.g., eq. [2])- Using direct numeri- 
cal simulations, we have shown that live dark matter 
satellites with realistic internal structure and mass ra- 
tios in the range relevant for galaxy formation (0.02 < 
Mgat/Mhost < 0.3) take much longer to merge than 
the corresponding point-mass satellite results, which are 
typically used in semi-analytic galaxy formation mod- 
els (e.g., Figs. |3][S]) . This difference is primarily because 
the satellites undergo significant mass loss as they sink 
deeper in the potential well of the host halo. Including 
a stellar bulge in both the satellite and host changes the 
merger time by < 10% for typical orbital circularities (see 
Fig. [5]). We find that the surprisingly simple fitting for- 
mula given in equations ^ and ^ provides an accurate 
(< 10% error) fit to the numerically determined galaxy 
merging times over the range of mass ratios, orbital cir- 
cularities, and orbital energies we have considered. The 



dependence of T mcrg0 on mass ratio, energy, and angu- 
lar momentum implied by this fit is all somewhat dif- 
ferent from standard assumptions in the literature (see 
£|3.3[) . We have chosen to test our results for a range of 
parameters that are cosmologically interesting, e.g. for 
satellite-to-host mass ratios for which the satellite under- 
goes significant orbital evolution within a Hubble time, 
and for orbital parameters that are relvant according to 
cosmological dark matter simulations. 

The results of this paper should be relevant to a wide 
range of problems in galaxy formation and evolution. For 
example, halo-halo merger rates as a function of mass 
ratio Mgat/Mhost from extended Press-Schechter the- 
ory or cosmological simulations show that a significant 
amount of dark matter mass is accumulated via accre- 
tion of lower-mass ha los with 0.02 < M„^JM^ nRt < 0.3 
(e.g., eq. [5] & Fig. IfJl lLacev fc ColdfloM IZentnerlfeoOTt 
Fakhouri & Ma 2007, in prep). Understanding the evo- 
lution of the galaxy population in such halos thus re- 
quires an accurate understanding of the effects of dy- 
namical friction in precisely the mass range where the 
point mass approximation is inadequate. To give a few 
concrete examples, we have shown that standard merg- 
ing timescales in the literature overestimate the growth 
of stellar mass by satellite accretion by ~ 40%, with the 
extra mass gained in low mass ratio mergers (Fig. |6]) . In 
addition, we have quantified the tendency for satellites 
that accrete onto a central galaxy to have lower angular 
momentum than average: the peak in the circularity dis- 
tribution shifts from i=s 0.5 to « 0.42 when one considers 
only satellites that merge onto a central galaxy within a 
Hubble time (Fig. [7j). 

Several limitations of our calculations should be noted. 
For orbits with very low circularity (77 < 0.2), our results 
are not applicable because such orbits pass sufficiently 
close to the center of the host halo on their first pericen- 
tric passage that interaction with the central galaxy is 
important and will decrease the merging timescale. All 
but two of our 28 simulations neglect the effect of the 
central galaxy in a given dark matter halo on the merg- 
ing timescale. As Fig. [2] shows, this is a reasonable ap- 
proximation when considering merging from the virial 
radius (as we have done), but would be less applicable 
for studying the detailed evolution of satellites at small 
radii in a dark matter halo (e.g. lOstriker fc Hausmanl 
[1971 . Finally, our results do not take into account drag 
on ambient gas, which m ay be importan t in cluster and 
group environments (e.g.. IOstrikeil[l999l ). 
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